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We compute the hydrodynamic torque on a dumbbell (two spheres linked by a massless 
rigid rod) settling in a quiescent fluid at small but finite Reynolds number. The spheres 
have the same mass densities but different sizes. When the sizes are quite different the 
dumbbell settles vertically, aligned with the direction of gravity, the largest sphere first. 
But when the size difference is sufficiently small then its steady-state angle is determined 
by a competition between the size difference and the Reynolds number. When the sizes 
of the spheres are exactly equal then fluid inertia causes the dumbbell to settle in a 
horizontal orientation. 


1. Introduction 

This paper is concerned with the question of how fluid inertia affects the settling of a 
particle under gravity in a fluid at rest. Understanding how particles settle in a fluid is 
of fundamental importance in many different scientific problems. An important example 
is sedimentation in suspensions of non-spherical particles. Homogeneous suspensions of 


settling spheroids are unstable to the formation of spatial patterns (Koch & Shaqfeh 1989 
Mackaplow fc Shaqfeh|1998[ Tornberg fc Gustavsson|2006[[Metzger et a/.|2007| Shin et al. 


DahlklidpOll Zhang et aZ.||2014|. The question is reviewed by Guazzelli & Hinch 


2009 


(2011). Other examples are plankton dynamics in the unsteady ocean (Jonsson 1989 


Guasto et al. 2012 Gustavsson et al. 2016), particle deposition (Charm et al. 2004), 


and the settling of ice particles in clouds, a process relevant for the onset and nature of 
precipitation ( Pruppacher fc Klett||2010 ). 

When the particle is very small, the particle dynamics can be computed in the creeping- 
flow limit, assuming that inertial effects are negligible. For larger particles, when the 
particle Reynolds number is finite, then fluid inertia must be considered. It is known, 
for example, that weak fluid inertia causes a sphere to settle more slowly. This is a 
consequence of the disturbance flow created by the sphere at small particle Reynolds 
number. The corresponding disturbance problem (the ‘Oseen problem’) was solved by 
|Proudman fc Pearson (19571 using matched asymptotic expansions. Fluid inertia also 


affects the transient settling of a sphere, before the steady state is reached ( |Sano||1981 
Lovalenti fc Brad^|1993 Mei||1992 ). 


For non-spherical particles less is known about the effects of weak fluid inertia. Only 
quite recently it was determined from first principles how weak fluid inertia modifies 
the torque on a spheroid in a simple shear (Subramanian fc Koch|[2005 Einarsson et al. 


20156|a Candelier et fl/.|[2015 Rosen et al. 2015). 

How does a non-spherical particle settle in a fluid at rest? Consider a small needle 
released at an inclination angle a with the horizontal. If fluid inertia is negligible then 
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Figure 1. Illustration of the dumbbell and the notation used in this article. The centres of 
two spheres (A and B) of diameters b + 5b and b are linked by a massless rigid rod of length 
a. The centre-of-mass of the dumbbell is denoted by C, and its velocity by v. The angle of the 
centre-of-mass velocity with the direction of gravity (negative e 2 -direction) is denoted by jd. As 
drawn the angle is negative, /3 < 0. Since we consider a steady state where the angular velocity 
of the dumbbell vanishes (see text), the centre-of-mass velocities of both spheres are also v. The 
angle of inclination of the dumbbell (defined by the unit vector n) with respect to the plane 
orthogonal to gravity (the ei-ea-plane) is denoted by a. 


the needle continues to settle at this angle if it is fore-aft symmetric. Its centre-of-mass 
moves at an angle /3 with the direction of gravity (Fig.[^. This angle is given by 


tan/3 = 


sin(2a) 
cos(2a) — 3 


( 1 . 1 ) 


Eq. (1.1) is obtained in the creeping-flow limit for a slender body with fore-aft symmetry, 


as discussed by Happel & Brenner (1983). See Eq. (5-11.33) in this book. 


Khayat & Cox (1989) analysed how weak fluid inertia affects the settling needle in 


the slender-body approximation. They found that a needle with fore-aft symmetry ex¬ 
hibits two equilibrium orientations: vertical (a = 7r/2) and horizontal (a = 0). Only the 


horizontal equilibrium is stable. Khayat & Cox (1989) concluded that the needle turns 
as it falls until it reaches horizontal orientation, a = 0, and then it continues to fall 
horizontally. This implies that its centre-of-mass velocity is aligned with the direction of 
gravity, that is /3 = 0. 


Dabade et al. (2015) computed how weak fluid inertia affects the orientation at which 


spheroids with arbitrary aspect ratios settle in a quiescent fluid. They found that when 
settling in a fluid at rest, weak fluid inertia causes an oblate spheroid to align its axis of 
rotational symmetry with the direction of gravity. 

Now consider a slender body that is rotationally symmetric but does not possess fore- 
aft symmetry. In the creeping-flow limit one expects the body to turn as it settles until it 
is oriented vertically (a = 7r/2). But fluid inertia exerts a torque that tends to turn the 


body towards the horizontal. As mentioned by Khayat & Cox (19891 this may cause the 


slender body to fall at an equilibrium angle that depends upon the degree of asymmetry of 
the body. In this paper we calculate this angle for a dumbbell without fore-aft symmetry. 


We use the method of reflection (Kim & Karrila 1991) to calculate the torque on the 
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dumbbell, taking into account convective fluid inertia to linear order in the Reynolds 
number Re;, = bv jv. Here b is the diameter of the needle, v is its centre-of-mass velocity, 
and v is the kinematic viscosity of the fluid. 

To conclude this Introduction we briefly comment on the wider context of this work. 
More generally the settling of particles in turbulent flows is of interest. Maxey (1987) 
showed that turbulence can enhance the settling speed of spherical particles. Recent 
experimental and numerical results on the settling of spherical particles in turbulence 
are discussed by Good et al. (2014). The dynamics of neutrally buoyant non-spherical 
particles in turbulence has been investigated by a number of authors (Parsa &: Voth|20T4 


Gustavsson et aL|2014 Byron et al. 2015 Soldati & Voth 2017), but little is known about 


the settling of heavy non-spherical particles in turbulence. 


2. Formulation of the problem 

We consider a simple model system: a dumbbell composed of two spheres, A and B. 
Their centres are linked by a massless rigid rod of length a (Fig. [^. Both spheres have 
the same mass density Pp so that the centre-of-mass of the dumbbell coincides with its 
geometrical centre. This means that the dumbbell experiences no gravitational torque 
with respect to its geometrical centre. But we mention that the effect of gravitational 
torques on particles with inhomogeneous mass densities (‘gyrotaxis’) is important for the 
motion of certain plankton species ( Durham et a7|[2013 Zhan et '^|2014 ). 

Sphere A has diameter b -|- Sb. Sphere B has diameter b. We take 6b > 0 and define two 
dimensionless parameters that characterise the shape of the dumbbell 


0 


and n = - . 

b 


( 2 . 1 ) 


Here A measures of the asymmetry of the dumbbell and k defines its aspect ratio. We 


note that Khayat & Cox (1989) define their parameter k differently. In the following it 


is assumed that the dumbbell is slender and that its asymmetry is small: 


K 3> 1 and A 1 . 


( 2 . 2 ) 


We do not consider the transient dynamics of the dumbbell but analyse its steady-state 
dynamics where the angular velocity of the dumbbell vanishes: 

u) = 0 . (2.3) 


In this case the two spheres settle with the same constant velocity v. We write: 

u = ?;sin/3 ei — ucos/3 02 . (2.4) 

Here jS is the angle between the velocity vector and the direction of gravity (the negative 
e 2 -direction. Fig. [^, and v is the magnitude of this velocity. To completely determine 
the steady-state dynamics of the dumbbell it remains to specify its angle of inclination 
a with the plane orthogonal to gravity (Fig. [^. We write: 

fi = cos{a) 01 -I- sin(a) 02 . 

The three unknowns of the problem {v, a and /3) are determined by three equations that 
are derived from the steady-state condition, namely that the force on the dumbbell and 
the torque with respect to its centre-of-mass C must vanish in the steady state. The force 
condition reads: 


0 = /a + /b 


( 2 . 5 ) 
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Here and /g are the forces acting on spheres A and B, and xa and xb are the 
position vectors of their centres. The torque condition 


0 = ta + (a^A - xc) X /a + tb + (a^B - a;c) x /j 


( 2 . 6 ) 


is a sum of the hydrodynamic torques centered on the spheres, ta and Tb, and the 


torques with respect to the centre-of-mass xc of the dumbbell. Eqs. (2.5) and (2.6) 


provide three independent conditions because the forces lie in the plane spanned by ei 
and 02 . The first equation gives two conditions, the second just one. It turns out that 


the torques ta and tb are negligible (Section 3). In this case it follows from Eqs. (2.5) 


and (2.6) that the forces /a and /b must align with the vector a:A — a;B. 

The spheres are subject to gravity and to the forces exerted by the fluid. The main 
difficulty lies in determining the force exerted by the fluid on each sphere. To this end 
we must solve the Navier-Stokes equations governing the fluid velocity and pressure. We 
note that the partial time derivative of the fluid velocity evaluates to zero in the frame 
translating with the centre-of-mass of the particle, since the dumbbell is assumed to settle 
with a time-independent velocity. In this frame of reference the Navier-Stokes equations 
read: 


V • 16 = 0 , p{{u ■ V)i6 = — Vp -I- pAu + pf g . 
The boundary conditions read: 


0 for X G = S^A U oS^B 


and 


u —> —V 


as 


oo. 


(2.7) 


( 2 . 8 ) 


Here u is the fluid velocity as seen by the particle, p = piV is the dynamic viscosity of the 
fluid, V is the kinematic viscosity, pf is its density, p is pressure and g is the gravitational 
acceleration pointing in the —e 2 -direction, and J^a and 5^b are the surfaces of the spheres 


A and B. We follow Maxey & Riley (19831 and decompose the fluid velocity and pressure 
as follows: 

u =—V + and P = {pt g ■ x + Pao) + ■ (2.9) 

Here —v is the undisturbed fluid velocity in the frame translating with the centre-of-mass 
of the particle, pi g ■ x + p^o is the (undisturbed) hydrostatic pressure, and and p^ 
are corrections to the fluid velocity and pressure due to the disturbance caused by the 
dumbbell. 

The forces /a and /b are given by 


/a = (toa - TOfA)^ + / and /b = (wb - Wfe) 9 + / crMs . 


( 2 . 10 ) 


The terms proportional to g are forces due to the undisturbed pressure (Archimedean 
forces) and gravitational forces. Here tua and tub are the masses of the spheres and rufA 
and mfB are the equivalent fluid masses. The remaining terms in Eq. ( 2.10[ ) are forces due 
to the disturbance. They are integrals over the fluid-stress tensor cr^ = —pH-f2/iS^ of the 
disturbance. Here is the symmetric part of the matrix of gradients of the disturbance- 
flow velocity u^. The integrals are over the surfaces and S^b of the spheres, ds is the 
surface element defined by the outward unit normal. The torques r a and tb are given 

by 


ta = 


[ (a; - cca) X o-Ms , 

Jy’A 


( 2 . 11 ) 


and Tb = / (a; — ojb) x cr ds . 

/jXa Jy-B 

To determine cr^ the disturbance problem must be solved. In the steady case the distur¬ 
bance velocity must satisfy: 


V • 16^ = 0 , Pf [(t6^ — d) • V] 16^ = — 


Vp^ -b pAu^ , 


( 2 . 12 ) 
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with boundary conditions 

= V for X G and 0 as |a;| —>■ oo . 


Eqs. (2.5), (2.6), and (2.10) to (2.13) constitute the problem to be solved. 


(2.13) 


3. Creeping-flow limit— method of reflection 

In order to determine the forces due to the disturbance flow we must determine cr^ by 


solving Eq. (2.12) with the boundary conditions (2.13). 


Consider how the orders of magnitude of the convective terms on the left-hand side 


of Eq. (2.12) compare to the magnitude of the viscous terms on the right-hand side of 


the same equation, in the vicinity of each sphere. It follows from the boundary condition 
Eq. (2.13) that the disturbance velocity is of order of v. The length scale is given by 


b. This allows us to estimate: 


Pf{u^ ■ V)m^ P{{v ■ V)m^ 


p,Au^ 


pAu^ 


O(Reb) . 


(3.1) 


In the creeping-flow limit. Ret = 0, these terms are negligible so that Eq. (2.12) becomes 


the steady Stokes equation determining the hydrodynamic interactions between the two 


spheres. We compute these interactions using the method of reflection (Kim & Karrila 
[1991 ). The method of reflection is applied iteratively. First the sphere A is considered as 
if this sphere was alone in a fluid at rest. The disturbance flow at x can then be written 


as: 


u 


1 ( 1 ) 


®) = - 1 + 


b^{i + xy 

24 


A 1 G(a; - xa) ■ 


(3.2) 


(the superscript in parentheses denotes the number of the iteration). The first term is 
the elementary solution of Stokes equation with a point force (Stokeslet). The second 
term is a short-range correction due to the finite size of the particle. Moreover 


/a^^ = —Z-npbil -I- A) i) 


(3.3) 


is the Stokes drag experienced by sphere A, xa is its position, and G(y) is the Green 
tensor of the Stokes equation: 




Stt/j, V 


r = |y| 


(3.4) 


The minus sign in Eq. (3.2) arises from the fact that the sphere is seen by the fluid as a 


point force, with intensity equal to the Stokes drag experienced by the sphere, but with 


the opposite sign. The terms in Eq. (3.2) that account for the finite size of the particle 


are of order (I-|-A)^6^'(;/r^. These terms thus contribute to the force acting on the sphere 
B at order at higher order than considered in our calculation. In the following we 
neglect these terms. 

Second the sphere B is introduced into the flow field Ma^\ disregarding its effect on 
sphere A. To determine the force acting on sphere B we use the first Faxen law (see for 
Happel fc Brenner||1983[) which leads us to 


instance 


- -?>T^pb[v - 


1 ( 1 ) 


62 


(®b) - 


1 ( 1 ) 


(aiB)] • 


(3.5) 


The term involving the Laplacian is neglected in the following since it is of order k 
Third, consider the torque Tb. In the creeping-flow limit the torque centred on sphere 
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B is given by 

~ Trfib^ftf{xB) where fl{ = x . (3.6) 

Since J7f(a;B) ~ bv/a"^ we see that this torque is smaller than the term (a^e — a:c) x 
by a factor of k~^. 

Fourth, in the same way force and torque acting upon sphere A are computed. Taking 
the results together we find to order k~^: 


/a = (wa — mf\)g — 37 r/r&(l + A)(i; — Sirfib G(—a) • v) + 0 {k , 

/b = (tob — rn{B)g — SiTf^b^v — 37 r^ 6 (l + A) G(a) • v) + 0 {k~^) . 


(3.7) 


Here a = xb — xa- Inserting these expressions into Eqs. (2.5) and (2.6) and noting that 
the torques ta and tb are negligible we can determine a, v, and /3 in the creeping-flow 
limit. 


3.1. Symmetric dumbbell (X = 0) 


For a symmetric dumbbell the torque condition (2.6) is always satisfied in the creeping- 
flow limit. It follows that all inclination angles a are steady-state solutions. This well- 
known result ( Happel fc BrenrLei]|1983 ) means that the dumbbell continues to fall at its 
initial inclination. For a symmetric dumbbell the force conditions (2.5) are readily solved 


to determine the centre-of-mass speed v and the angle /3 as functions of a. We find: 


tan /3 = 


sin(2Q;) 


cos( 2 c 


o 16 

+ 3-y» 


-Fs^ + O(0. 


(3.8) 


We have expanded Eq. (13.81) for large values of k because our perturbative solution 


assumes that k, is large. In this limit the angular dependence in Eq. (3.8) is qualitatively 
similar to that in Eq. (Hi, derived for a thin needle with fore-aft symmetry ( |Happel~fc 
Brenner 1983). According to Eq. ( |3.8| ) the angle /3 is negative for 0 < a < 7 r /2 and it 
tends to zero as a —>■ 0 and a —>■ 7 r/ 2 . An important difference between Eqs. (1.1) and 


(3.8) is, however, that the angle (3 is independent of the aspect ratio in the case of the 


needle, while it does depend on k for the dumbbell. According to Eq. (3.8) the angle /? 
tends to zero in the limit k —>■ oo. This is simply a consequence of the fact that there are 
no hydrodynamic interactions between the spheres in this limit; they must fall as if they 
were independent, each at the same terminal velocity. 

For the centre-of-mass speed we find: 


V = 


mfB(7 - l)g 

Unbyt 


1 -f 


6 — 3cos(q:)" 
8k 


+ 0{k-^) 


(3.9) 


where 7 = uiBlrrifB- From (3.9) we infer that vertical orientation [a = tt/2) yields the 
fastest terminal velocity. We also see and that, in general, the dumbbell settles faster 
than either sphere A or B alone. This is due to the disturbance flow. 

3.2. Asymmetric dumbbell (X^ 0) 

Now consider an asymmetric particle. When A is not zero then Eqs. (2.5) and (2.6) admit 
the solutions: 


and 


a = ±tt/2 , /3 = 0 

^£3(7-1)5/- 3 


SiTfib 


(i + £ + a ) +0{X^,k-^,Xk-^) . 


(3.10) 

(3.11) 
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The asymmetric dumbbell settles vertically in the steady state. From the torque equation 


(2.6) it follows that a = — 7 r /2 (the smaller sphere first) is unstable. By contrast, a = 7 r /2 


is stable. So when an asymmetric dumbbell is released in fluid at rest then it tends to 
settle vertically, the larger sphere first. 


4. Effect of fluid inertia 

In this Section we incorporate weak effects of fluid-inertia on the dynamics of the 
dumbbell settling in a quiescent fluid, assuming that Re;, small but finite. As mentioned 


in the Introduction, Khayat & Cox (1989) investigated the effect of fluid inertia on the 


settling of needle with fore-aft symmetry (see also Dabade et al. 2015). What happens 


when fore-aft symmetry is broken? To answer this question we need to determine the 
disturbance flow produced by a single sphere when convective inertia is taken into ac¬ 
count. The disturbance flow produced by a single sphere corresponds, to leading order 
and in the vicinity of the sphere, to a Stokes flow plus a small correction. Since the 
Stokes solution scales as vb/r it follows that far from the spheres the convective term 
P[{u^ ■ V)m^ is always smaller than the viscous term pAu^ and can be neglected. But the 
convective term pf{v ■ V)m^ decays more slowly than the viscous terms and must balance 
the viscous terms at a distance 6/Ref,, the Oseen length. At this distance this convective 
term can no longer be neglected. In this way we are led to the Oseen equations. 

4.1. Elementary solution of the Oseen equations 
We use an iterative procedure, as in the previous section. The first step is to find an 
elementary closed-form solution of the Oseen equations. Keeping in mind that the two 
spheres are assumed to be far from each other (k ^ 1 ) the rapidly decreasing terms 
with respect to r in the solution of the Oseen equations can be neglected. We therefore 
consider the equations 


V • M = 0, —pf{v ■ V)m = — Vp + pAu + f6 . 


(4.1) 


where the particle is represented by a point force {S{x) stands for the Dirac delta func¬ 
tion). The solution of Eqs. (4.1) can be found in the literature (Pignatel et al. 2011), 
see also (Happel & Brenner (1983), p. 79). It can be obtained as follows. Applying the 


Fourier transform 


we find 


u = Fu = / u{y) exp(—z k ■ y) dy 


k ■ u = 0, —ipf {v ■ k)u = —ikp—pk^u + f. 


(4.2) 


(4.3) 


As usual, the pressure is determined by projecting equation (4.3) along fc. This yields 
—ip = —f • and performing a partial fraction decomposition we are led to 

/ f f k 


u = 


i}i[ — _ ^ _ 

kf {pk — ipiv) ■ k 


(pk — ipfv) ■ k \pfV ■ k 
This result simplifies: since the vectors f and v are collinear we have 

\PiV ■ kj pfv ’ 

where / = |/|. Now we apply the inverse Fourier transform. Using the relations 


= 


1 


47r r 


and F 


1 


{pk — ipfv) ■ k 


dirp r 


(4.4) 


(4.5) 


(4.6) 
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we find 


6 xp[- + ^ ■ y)] J, f ^ l-exp[ - ^{vr + v-y) 

ATTfj, r pfV 47r r 


Expanding the derivatives in Eq. (4.71 yields 


_exp[- ^{vr + v-y)] ^ 

u — _ j 


+ 


Sttp r 


(4.7) 


(4.8) 


- + ■ y) } 


y 


X _ 

PfV 47r 


The next step is to apply the method of reflection, with Eq. (4.8| instead of Eq. (3.2). 


4.2. Method of reflections 

As a first step the flow produced by the sphere A is considered as if it were alone in 


a fluid at rest. The result is given by Eq. (4.8) in which / is replaced by 37r5(l + A)t; (that 
is the Stokes drag, with a minus sign). The sphere B is then introduced into this flow. We 
note first that the flow produced by the sphere A is non-uniform. In the previous section 
we discussed the fact that Faxen corrections are negligible in our problem. One may ask 
whether the non-uniformity of may give rise to a Saffman lift force acting on the 
sphere B. It turns out that, in the equations governing the flow produced by the sphere B, 
effects of the convective terms due to the non-uniformity of are negligible compared 


u 


1 ( 2 ) 


to those arising from the Oseen convective terms of the form pf { [d —(^^b)] • V } 

The lift force can be therefore be neglected and the equations which govern the fluid 


velocity take the same form as Eqs. (4.1) but where the velocity v has to be replaced by 


the relative velocity v — M^^^^(a:B)- The force acting on the sphere B (see for example 
Proudman fc Pearson|[l95^ is then given by 


= -37rp6(l -I- [t> - ^^^^(a^B)] + O (k AReb, Re;,K 


(4.9) 


(note that depends on Ret,). Similarly, the disturbance force acting on the sphere 

A reads as 


= -37r/i6(l -k A)(l -k ^Ret,) [v - ■ 




(oja)] + O (k a Ret,, Ret,K ^) . (4.10) 


4.3. Results 


Inserting the disturbance forces (4.9), (4.10) into Eqs. (2.10) and solving Eqs. (2.5) 


and (2.6) yields the desired solution. In general this must be done numerically. Here we 


consider the limit tcRet, ^ 1. In this case we can obtain explicit expressions by expanding 


the Oseen solution (4.8). In this way we find that the angle of inclination tends to the 
equilibrium value 




arcsm | 

TT 
2 


3 Ret, 


, 16 A 

when —-— 
3 Ret, 

otherwise. 


< 1 , 


(4.11) 


As expected it follows from Eq. (4.11) that a sufficiently asymmetric dumbbell falls in 


vertical orientation. A perfectly symmetric dumbbell, by contrast, settles in the horizon¬ 
tal orientation, like a fore-aft symmetric slender body (Khayat & Cox [1989 ). This is a 
consequence of the fact that the Oseen solution breaks the symmetry of the Stokes so¬ 
lution. But when the effect of breaking fore-aft symmetry and of convective fluid inertia 
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balance then the dumbbell falls at the equilibrium inclination determined by Eq. (4.11). 
For Rch = 0.1, for example a slight asymmetry of one percent, A = 0.01, causes the 
dumbbell to settle at an inclination angle different from zero. The vector n forms an 


angle of roughly 30 degrees with the horizontal according to Eq. (4.11), as schematically 


shown in Fig.[^ One expects slender particles to exhibit a qualitatively similar behaviour 
and this could be tested experimentally by observing a small rod settling in a quiescent 
fluid at small Reynolds number. 


Eq. (4.11) shows that the angle is independent of k when nReb small. This is a conse¬ 


quence of a property of the Oseen solution that is particular to the dumbbell. Expanding 
the Oseen solution assuming that a is of the order of or smaller than 5/Reb, we deduce 
that there is a uniform contribution to the fluid velocity that is proportional to Ret,. This 
means that the inertial torque is proportional to a. Since the contribution to the torque 
from the particle asymmetry is also proportional to a it follows that a is independent of 
K. A slender body behaves slightly differently, since the torque has a different dependence 
upon a. See Eqs. (6.12) and (6.22) in ( Khayat fc Cox|[l9^ ). When /tRe{, <C 1 the angu¬ 
lar dependence of the inertial torque is the same for the dumbbell as for a thin needle, 
proportional to sin( 2 a). 

At the order considered here, the angle /3 of the centre-of-mass velocity is found to 


remain unchanged compared with the creeping-flow limit (3.8). The settling speed is 
given by 


mmh - l)g 

SirbyL 


1-b 


6 — 3cos(a)" 
8k 


-f A — — Re^ 


+ 0(e^). 


(4.12) 


Here denotes any quadratic combination of the small parameters k~^, A, or Re;,. 
Eq. (4.12) is consistent with the results obtained in the creeping-flow limit in Sect ion 
In the limit Reb —>■ 0 Eq. (4.12) reduces to Eq. (3.9) for a ^ 0, A = 0, and to Eq. (3.11) 
for a = 7 r/ 2 , A 7 ^ 0. 


5. Conclusions 

We computed the hydrodynamic torque on an asymmetric dumbbell settling in a qui¬ 
escent fluid, assuming that the Reynolds number Ref, is small but finite. The two spheres 
have the same mass densities but different sizes. This asymmetry gives rise to an addi¬ 
tional contribution to the torque that may balance the contribution of convective fluid 
inertia. In this case the dumbbell settles at an equilibrium angle that is determined by 
this balance. This prediction is expected to qualitatively hold more generally for particles 
of other shapes with broken fore-aft symmetry. When bcRet is small the equilibrium angle 
of the settling dumbbell is independent of the value of the aspect ratio k. This property 
is particular to the dumbbell, it does not hold for more general bodies. In an indepen¬ 


dent study, Roy et al. (2016) have analysed the settling of small rod-like and ramified 


particles without fore-aft symmetry that settle in a quiescent fluid, using slender-body 
approximations and experiments. It will be of interest for future work to determine the 
full range of particle shapes and flow parameters at which fore-aft asymmetry produces 
an equilibrium sedimentation angle that is neither horizontal nor vertical. It is of par¬ 
ticular interest to answer the question: how robust are results derived assuming fore-aft 
symmetry? 

Our results were obtained using the method of reflection. This method generates a 
perturbation expansion in the inverse of the aspect ratio k of the dumbbell. The method 
makes it possible to systematically treat the dynamics of more general assemblies of 
spheres, partially linked. It allows, for instance, to compute the dynamics of two dumb- 
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bells settling together, affecting each other but not coming too close to each other. The 
method used here is not restricted to quiescent flows. Applying the method of reflection 
for a neutrally buoyant symmetric dumbbell in a simple shear gives Jeffery’s equation 
( Jeffery|1922 ), and it may be possible to treat the effect of fluid inertia upon the angular 
dynamics of the dumbbell using the method of Section and the reciprocal theorem 
(Subramanian fc Koch|[2005 Einarsson et fl/.|[2015a|6 ). 
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